Semiclassical structure of chaotic resonance eigenfunctions 
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We study the resonance (or Gamow) eigenstates of open chaotic systems in the semiclassical 
limit, distinguishing between left and right eigenstates of the non-unitary quantum propagator, and 
also between short-lived and long-lived states. The long-lived left (right) eigenstates are shown to 
concentrate as h — > on the forward (backward) trapped set of the classical dynamics. The limit of 
a sequence of eigenstates {ip(h)}h^o is found to exhibit a remarkably rich structure in phase space 
that depends on the corresponding limiting decay rate. These results are illustrated for the open 
baker map, for which the probability density in position space is observed to have self-similarity 
properties. 
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In closed systems the two most fundamental semiclas- 
sical properties are that the mean density of states is 
given by the Weyl law |l| , which associates to each quan- 
tum state a Planck cell in the available region of phase 
space, and that in classically chaotic systems the station- 
ary states have Wigner functions which are semiclassi- 
cally uniform over the energy shell 2] , in agreement with 
the quantum ergodicity theorem 3]. It is remarkable 
that it is still not known in general how these fundamen- 
tal properties extend to open (scattering) systems. 

In open systems the lack of unitarity of the quantum 
propagator gives rise to non-orthogonal decaying eigen- 
states with complex energies (resonances), the imaginary 
parts of which are interpreted as decay rates. In the case 
of open chaotic systems the classical mechanics is struc- 
tured in phase space around fractal sets associated with 
trajectories that remain trapped for infinite times, ei- 
ther in the future (forward-trapped set A + ) or in the 
past (backward-trapped set A"_). The mean density of 
resonances is believed (but not in general proved) to be 
determined by the fractal dimension of the invariant set 
Kq = K + n A_,the classical repeller. This is the frac- 
tal Weyl law 0,0,01. (Note that this is different to the 
resonance statistics in weakly open systems, for which 
the size of the opening vanishes in the semiclassical limit 
0) 

Much less is known about the resonance (or Gamow) 
eigenstates. These are important in many areas of 
physics H and chemistry 0, because they have marked 
influence on observable quantities such as scattering 
cross sections and reaction rates (they are a compo- 
nent of the Siegert pseudostates basis in terms of which 
the scattering wavefunctions and S matrix, for exam- 
ple, can conveniently be expanded [Sj)- Following the 
well established idea that in the semiclassical limit time- 
independent quantum properties should be related to 
time-independent classical sets, it is natural to expect 
that long-lived eigenstates of open systems should be de- 
termined by the structure of K + and A"_ . This was tested 
numerically for some right eigenstates of the open kicked 



rotator in [lfj , where the term 'quantum fractal eigen- 
states' was coined. 

We here significantly extend the notion of quantum 
fractal eigenstates in several new directions. First, we 
draw the important distinctions between left and right 
eigenstates of the non-unitary propagator, and between 
states that are 'short-lived' and 'long-lived' with respect 
to the Ehrenfest time. Second, we show that in the semi- 
classical limit the long-lived left eigenstates concentrate 
on A + , while the long-lived right eigenstates concentrate 
on A'_. In chaotic systems the eigenstates thus inherit 
the intricate fractal structure of the underlying classical 
trapped sets (this property has also been observed by 
Nonnenmacher and Rubin ^3). Third, we find that in 
the semiclassical limit the eigenstates have a rich struc- 
ture that reflects the self-similarity of the sets K± and 
depends explicitly on the limiting decay rate. A semiclas- 
sical formula is derived for the weights of the quantum 
eigenstates in different regions of phase space. This for- 
mula is consistent with the numerical results of ^(J m the 
special case when the decay rate is equal to the classical 
escape rate. 

We illustrate our results for an open quantum baker 
map. In this system we have, in addition, that in the 
semiclassical limit long-lived right eigenstates have frac- 
tal support in momentum space and are self-similar in 
position space (and vice- versa for left eigenstates). Fi- 
nally, we discuss a system, the Walsh-quantized baker 
map, for which exact results corroborate our arguments. 

Many recent studies of open quantum systems have fo- 
cused on maps of the torus T, viewed as a phase space, 
which is opened by removin g a strip parallel to the mo- 
mentum direction 0, 0, llH. Il4l fl5j | . A physical motiva- 
tion is the "bounce map" defined for billiards, in which 
case the opening corresponds to attaching a perfect lead 
to the billiard. The quantization of the classical map IA 
is a unitary matrix U, acting on a Hilbert space of di- 
mension N. For torus maps this dimension plays the role 
of an effective Planck's constant, h — (2ttN)~ 1 , and the 
semiclassical limit thus corresponds to N — > oo. If the 
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opening O occupies a fraction M/N of the total area in 
phase space then quantum mechanically the open map 
corresponds to a non-unitary matrix U = UIL, where II 
is a projector onto the complement of the opening. The 
result of multiplying by II is to set M columns of U equal 
to zero. Since the matrix U is not unitary, we must dis- 
tinguish between its left and right eigenstates, 

U\^)=z n \^), (^\U = z n (^\. (1) 

We will assume = (*^|*^) = 1. The eigen- 

values z n lie inside the unit circle in the complex plane, 
\z n \ 2 — e~ r " < 1 where T„ > is the decay rate. (Note 
that these are not scattering states, for which fractality 
has also been observed 0.) 

If O m = U m (0) denotes the mth image of the open- 
ing under the classical map, the forward-trapped and 
backward-trapped sets are defined, respectively, as 

oo oo 

K+ = T\\J 0_ m , K-=T\\J O m , (2) 

m— m— 1 

Let us also define the set of points which fall into the 
opening after m steps, but not earlier, 

TV£ = {x e 0- m , x £ C_ n , 0<n< m}. (3) 

We also define = {x 6 O m , x £ £>„, 1 < n < m} for 
m > 1 and TVl — 0\ . These sets are related by 

u- l (n™)\o = n™ + \ u{n m \o) = n m+ \ (4) 

with the convention that TZ , = O. 

The Ehrenfest time te plays an important role in quan- 
tum chaotic transport Till Il7j . essentially de- 
termining the quantum-to-classical crossover. This is 
the time it takes for a minimal wave packet to become 
larger than the opening, being therefore partially re- 
flected when leaving the system. For times shorter than 
te the quantum evolution of localized states can be ap- 
proximated by the classical evolution. For open chaotic 
maps te = A lnM, where M is the number of 'open 
channels' and A is the Lyapunov exponent UlEl. In 
the semiclassical limit quantum states that are essen- 
tially supported on some region A C 7£" of phase space 
will have a deterministic escape and will thus become 
eigenstates of U n+1 corresponding to eigenvalues that 
vanish as h — > 0. These are called 'short-lived' states 
[Til Il5| . The estimate that the fraction of such states 
is 1 — e~ TB / TD , which holds when the mean dwell time 
te = N/M is large, results in the fractal Weyl law [l5| . 
Short-lived states are intimately related to the short-time 
dynamics and are relatively insensitive to the classical 
trapped sets. In what follows we obtain several results 
for the semiclassical limit of long-lived eigenstates. 

One way to look at the eigenstates is to depict them 
in phase space using the Husimi representation, in which 



one associates with any eigenstate |\&|) (£ = L,R) the 
function H^(x) = | | 2 , where x = (q,p) is a point 

of phase space and \x) is a coherent state. From Q we 
have for the Husimi function of a right eigenstate 

\z n \ 2m \(x\^)\ 2 = \(x\U m \^)\ 2 . (5) 

For times much shorter than the Ehrenfest time the ac- 
tion of the quantum propagator on coherent states can 
be approximated by the classical evolution; in particular 
we have that as H — > 

ffi) m \x)RtO if xeTZ m . (6) 

Therefore, in the semiclassical limit becomes negli- 
gible on regions 1Z™ if \z n \ > is fixed and ra « tj. 
Since the Ehrenfest time grows as H — > 0, the function 
must concentrate on the backward-trapped set K- in this 
limit, 

H*(x)&0 if x^K- (\z n \ >0, ft-»0). (7) 

For a finite H there will be some leakage of |(^| V E'^')| 2 
into the regions O m , mostly for large m. This deviation 
should be more pronounced for states with small lifetime. 
The reasoning for long-lived left eigenstates is completely 
analogous. In the semiclassical limit they concentrate on 
the forward-trapped set K + . 

For finite h the Husimi functions are not supported 
on truly fractal sets, but as h — > finer classical struc- 
tures are revealed. Continuous time Hamiltonian sys- 
tems can be treated in much the same way. One just 
replaces \z n \ 2m by e~ TnT and the essence of the argu- 
ment remains. We thus expect Husimi functions that 
live on fractal sets in generic open chaotic systems in the 
semiclassical limit. 

Let us investigate the weight of the eigenstates on dif- 
ferent regions of phase space. Let ttq be the projector 
onto the opening. Our non-unitary propagator then sat- 
isfies U'U = 1 — ttq. Hence we find that the weight of the 
right eigenstates in the opening, 

«17T0|V]/«H1-|Z„| 2 , (8) 

increases with increasing decay rate. Let 7r m , m > 1, be 
the projector onto (this is an example of the general 
class of projectors introduced in (lRj). For short times 
the semiclassical approximation gives U^ir m U w tt to +i, 
and thus we obtain the remarkably simple relation 

(^\n m \^)^\z n \ 2m (l-\z n \ 2 ). (9) 

Notice that this vanishes for any fixed m as \z n \ — ► 1, 
thus if the quantum decay rate approaches zero as h — > 
the corresponding eigenstate becomes localized on the 
invariant set K$. 

It is interesting to note that the distribution of the 
eigenstates on phase space depends explicitly on the cor- 
responding decay rate, which we hold fixed as h — > 0. 



FIG. 1: (color online) Left panel: The average of the Husimi 
functions of the 100 longest- lived right eigenstates of the baker 
map, for N — 3 7 (intensity increases from blue to red). Right 
panel: the corresponding Wigner function average (in the 
white regions the function is non-positive). Note that the 
ft-scale is l/N « 0.005. 

This is an important difference with respect to the case 
of unitary evolution, when the semiclassical limit is the 
same for almost all sequences of eigenstates {tp(h)}h^o- 
We also note that the areas of the regions V,™ are pro- 
portional to e _m7c , where 7 r . is the classical escape rate. 
Therefore the prediction in [10( that a right eigenstate of 
an open quantum map with T n rs j c should be constant 
(up to quantum fluctuations) over the classical set K— is 
consistent with our expression. 

We now illustrate our results in a paradigmatic model 
of quantum chaos, the baker map, which is particularly 
convenient because of the relatively simple structure of 
its trapped sets. Its triadic form is defined as 

f(3g,f) if < q < i, 

U{q,p) = l{Zq-l,^) if !<<?<§, (10) 
[(3 9 -2,£±2) if §<<?<!. 

This system is uniformly hyperbolic, with Lyapounov ex- 
ponent A = In 3. The stable and unstable manifolds are 
parallel to the momentum and position axes, respectively. 
Its quantization is given by [19| 

U N = F N 1 di a ,g(F N/3 ,F N/3 ,F N/3 ), (11) 

where the subscript denotes the dimension of the ma- 
trix. Fn is a modified Fourier transform, (F]y) nm = 
-J=e~~^( n+1 / 2 ^ m+1 / 2 \ with the integers n,m running 
from to N — 1. The "1/2" factors are necessary in or- 
der to impose parity on the eigenstates. This form of the 
matrix Fn leads to antiperiodic boundary conditions. 

An open version of this system was introduced in 0, 
H^ . in the context of the fractal Weyl law, where the 
opening was taken as the middle vertical strip (sec also 
[2C| where a closely related map was studied) . To find the 
forward trapped set K + one may consider the backward 
propagation of this strip, and it is not difficult to see that 
the only points that remain are those belonging to Can x 
[0, 1), where 'Can' denotes the usual middle-third fractal 



FIG. 2: (color online) The weight (*^|7r m |*^) of right eigen- 
states of the baker map on the regions 1Z+ , for N = 3 6 . The 
solid line is the semiclassical approximation @. 

Cantor set. The backward-trapped set K- is [0, 1) x 
Can. Details of this construction can be found in 0, 
ll2T ]. The open quantum system is obtained by setting 
the middle third of the columns of C//v, corresponding to 
the strip, equal to zero. The kinematics is such that the 
position representation of right eigenstates is equal to the 
momentum representation of left eigenstates. 

In Fig.l we plot the average of H^(q,p) over the 100 
longest-lived states, for the case N = 3 7 . We see that it 
is reasonably concentrated on the backward-trapped set, 
although it is not able to resolve this set on the finest 
scales. More details of K- should be revealed for larger 
values of N. The right panel of Fig.l shows the averaged 
Wigner function, and we see that it resolves K- with 
considerably greater accuracy. 

In Fig. 2 we plot the weight of right eigenstates on the 
regions 1Z™, for N = 3 6 . The solid lines represent the 
semiclassical approximation 0. We see that the agree- 
ment is rather good, particularly for small eigenvalue 
modulus and small values of m. 

Let us now turn our attention to wave functions in 
position/momentum space. For the baker map, the pro- 
jection of K- onto momentum space is just the Cantor 
set, and the momentum representation of any long-lived 
right eigenstate must be supported on this fractal dust. 
In Fig. 3 we plot | ^ ^ (^?) | 2 , averaged over the 25 longest- 
lived states, for N = 3 7 . The fractal nature of the sup- 
port, within the resolution set by the uncertainty prin- 
ciple, is evident from the magnification by a factor of 3. 
In Fig. 4 we plot | (q\^ R ) | 2 averaged over a few right eigen- 
states, again for N — 3 7 . Panels a and b correspond to a 
relatively large decay rate, while c and d correspond to a 
smaller one. Within the available resolution both func- 
tions are remarkably self-similar, as can be appreciated 
from the magnifications. This observation is consistent 
with our theory (cf. @). 

Nonncnmachcr and Zworski have recently introduced a 
simplified version of the open quantum baker map, based 
on the 'Walsh-Fourier' transform, that can be solved ex- 
plicitly [l^. For this system several of the results pre- 
sented here are exact and can be rigorously proven. One 
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FIG. 3: Probability density in momentum space for right 
eigenstates, KpI^^)! 2 , averaged over the 20 longest-lived 
states (eigenvalue modulus ranging from 0.90 to 0.83). This 
function is approximately supported on the Cantor set, as can 
be seen from a magnification by a factor of 3 (right panel). 
TV = 3 7 and so the ft-scale is l/N w 0.005. 



particular feature is that the difference between short- 
lived and long-lived states is well defined. The former 
all have exactly null eigenvalue, and the phase space rep- 
resentation of the latter vanishes outside the finitely re- 
solved relevant trapped set. Equations © and © are ex- 
act for this system. We have also shown that self-similar 
properties analogous to those illustrated by Fig. 4 can be 
established rigorously and that -after either a spatial av- 
erage or a spectral average- the probability density of 
long-lived right eigenstates in the momentum represen- 
tation becomes constant on the Cantor set as N — > oo 
(for the standard quantization this is suggested by Fig. 3). 



A detailed account will be presented elsewhere 
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FIG. 4: Probability density in position space for right eigen- 
states, Kgl^^)! 2 , averaged over states with similar decay 
rates. Panels b and d are magnifications of a and c. The av- 
erage eigenvalue modulus is approximately 0.4 for a-b and 0.7 
for c-d. In both cases the self-similarity is striking. N = 3 7 . 



Summarizing, we have shown that in the semiclassi- 
cal limit resonance eigenstates of open chaotic systems 
concentrate on classical fractal trapped sets. We have 
derived a formula for their distribution on different ar- 
eas of phase space that depends explicitly on the quan- 
tum decay rate. For the baker map we have also found 
that the position and momentum representations of wave 
functions exhibit self-similarity. 
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